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In the proposed work, a defect cavity multi-layer Bragg reflector structure is 
proposed theoretically to find the presence of thyroid cancer cells in the 
given sample. The modelling, design and analysis of the sensor is performed 
using characteristic matrix method (CMM). Proposed structure has central 
defect cavity with 6 pairs of low and high refractive index layers on each 


side of the defect. To enhances the sensor sensitivity, the incident light in 


mid-infrared frequency range is used as input light source. The refractive 
index of normal and thyroid cancer cells is analysed for the performance of 
the sensor. The obtained Q factor and sensitivity of the sensor design is 3729 
and 2828 nm/RIU respectively. The proposed sensor is a best choice of 
optical sensor for the detection of thyroid cancer cells in the given test 
sample for accurate analysis in medical applications. 
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1. INTRODUCTION 

Thyroid, an important organ in the human body butterfly in shape and located below the Adams 
apple of the front neck region. It has two side lobes which is connected by isthmus [1]. Thyroid is rich in 
blood vessels and plays important role in voice quality. It generates various hormones which plays a 
significant role in the various parts of the body. Metabolism is one such process where the food digestion is 
converted into energy and is controlled by with the proper generation of thyroid hormones [2]. When thyroid 
generates too low or too high thyroid hormones, it causes hypothyroidism or hyperthyroidism respectively. 
The two specific and significant hormones generated by thyroid glands are thyroxine (T4) and 
triiodothyronine (T3) [3]. These two hormones are responsible for informing the body cells to properly utilize 
the energy available. In its normal working thyroid generates balanced level of these two hormones and 
control the body activity. The entire activity of thyroid is under the control of pituitary gland. It is located 
below the brain and in the center of the skull. This pituitary gland is responsible for sensing and controlling 
the amount of thyroid hormones in the blood stream. When pituitary gland senses the deficiency or increased 
amount of thyroid hormones, it automatically adjusts the hormones level by its own. This is done by 
generating thyroid stimulating hormone (TSH) hormones [4]. This TSH hormone instructs the thyroid gland 
to take some decisions to bring back the body with normal conditions. Thyroid disease is a general term used 
to mention the imbalanced level of thyroid hormones in the body-system. Increased thyroid hormones makes 
the body to consume more amount of energy too quickly. It also causes tiredness and reduction in the body 
weight. On counterpart, if thyroid generates too low level of hormones, the amount of energy consumed by 
the cells reduces and causes tiredness and makes to gain body weight. Thyroid disease may found in men, 
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women, teenagers, infants and elderly persons. It can be inherited and passed down through families. 
Hyperthyroidism can be found during the birth and may develop with the age. There are many reasons to 
have thyroid diseases viz, family history, anemia, type 1 diabetes, turner syndrome, with the age like more 
than 60. People with type 1 diabetes are more likely to get thyroid disease and compared to others as type 1 
diabetes is an autoimmune disorder [5]. Symptoms of hyperthyroidism includes, anxiety, weight loss, goiter, 
muscle weakness, nervousness, irritability [6] and that of hypothyroidism includes, gaining weight, fatigue, 
forgetfulness, hoarse voice, intolerance to cold temperatures [7]. For the year 2022 The American Cancer 
Society estimates about 43,800 cases of thyroid cancers with around 2,230 deaths [8]. There is about 0.6% 
increase per year in the death rate due to thyroid cancer. Studies concludes that women’s were 3 times more 
likely to attack for thyroid cancer than men’s. Blood test, imaging and physical examinations are the 
commonly used test to diagnose the thyroid cancer. These examinations have limitations like more sample 
requirement, less accuracy, more detection time, experienced lab technicians and costlier. Micro-electro 
mechanical systems (MEMS), micro-opto-electro mechanical systems (MOEMS) and optical sensors plays 
an important role in biomedical applications [9], [10]. Optical sensors are the best choice for sensing 
application as they require less sample, small in size and weight, more accurate, no electromagnetic 
interference (EMI) and can be fabricated using micro-machining techniques [11], [12]. 


2. THEORETICAL DESIGN AND METHODOLOGY 

Easy fabrication, simple design, quick and accurate sensing of a sample are the major factors that 
many researchers considered while designing sensors. The best choice to satisfy these criteria are 1D 
photonic crystals. By properly choosing the material used for designing the sensor, the properties of the light 
inside the PhC can be precisely controlled [12], [13]. Photonic band gap (PBG) is the unique property of PhC 
which indicates the frequency range prohibiting from the light propagation [14], [15]. By creating a defect in 
the structure, the light can be made to propagate in this frequency region. Structures containing stack of 
alternating layers of two different refractive index (RI) materials is called distributed Bragg reflector (DBR). 
It acts as 1D photonic crystal due to its periodic variation in the RI of the material in only one direction. 
As the light is incident on these structures, a part of the light is transmitted, and the remaining part is reflected. 
This can be observed in every boundary at the interface of two materials. Due to constructive interference, 
the reflected light combined and causes a very high reflection. Stopband is the frequency range where the light 
propagation is forbidden and hence exhibits the property of photonic bandgap [16]-[18]. Figure 1 shows the 
structure of 1D DBR consists of alternating layers of two different RI materials. The structure considered has 
N = 6 pairs of low (L) and high (H) RI material. 
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Figure 1. Proposed multi-layer DBR structure having N = 6 


To calculate the RI value of each material we can use the following formula: 


Ny = VE N, = VEL (1) 


Where, n; and ny are the RI of low and high layer respectively with the permittivity €; and €y. 

The incident electromagnetic wave experience multiple reflections as it propagates through DBR. 
For this multiple Bragg reflection, the optical length of each layer should be equal to the quarter wavelength 
of the incident light. 
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Where, l, and ly indicates the geometrical thickness of layer low and high RI material respectively 

Ao = = is the wavelength of incident light and c being the free space velocity of the electromagnetic wave. 
0 

Figure 2 shows the transmission spectrum of N = 6 pair DBR. It is clear from the transmission spectrum that, 


for a certain frequency range there is no transmission. This region is called photonic band gap (PBG) region. 


Transmission (%) 


0 ‘ 
10.5 11 11.5 12 12.5 13 
Wavelength (4m) 


Figure 2. N = 6 pair multi-layer DBR structure transmission spectrum 


In this work, the high RI layer (H) is made with germanium (Ge) material and the low RI layer (L) is 
made with zinc sulphide (ZnS) material with RI = 4 and RI = 2.19 respectively. For the devices operating in 
the mid-infrared wavelength range (3 um — 25 um), Ge and ZnS are the commonly chosen materials. They have 
the advantages of lower absorption coefficient, higher refractive index contrast and easy fabrication process. 
By using pumped CO2 lasers, these devices can be operated at THz frequency ranges [19], [20]. Sol-gel and 
spin coating methods make the DBR fabrication easy [21]—[23]. Reactive electron beam or chemical vapor 
deposition (CVD) methods of deposition can also be used for the fabrication of multi-layer structures. 
The imaginary part of RI of the materials used here is neglected as these materials have very low absorption 
coefficient at the frequency range used here. For the proposed design, mid-infrared wavelength of 10.6 um 
incident light is considered. Because of good molecular fingerprints in the mid-infrared region, it offers higher 
sensitivity of the sensor [24], [25]. The obtained geometrical thickness for low and high index layer is 1.21 um 
and 0.6625 um respectively. All the simulations were done using Matlab tool. 


3. CHARACTERISTIC MATRIX METHOD 

In characteristic matrix method (CMM) is used to analyse the interaction of incident 
electromagnetic waves with the layered media. The incident light is assumed to incident normally on the first 
layer and is transverse electric (TE) polarized. If the angle of incidence is denoted by 0 then, for TE wave, 
we have [13], [26]. 


du 


2 
— + (kĝn?cos?0)U = 0, TY + (k2n?cos?0)V = 0 (3) 


dz2 


The solution obtained for the (3) are: 


U(z) = Acos (konz cos@) + B sin (konz cos@) 
V(z) = + ficos 0 {B cos(kpnz cos@) — A sin (konz cos@)} 


The obtained characteristics matrix is: 


uo = cos (konz cos@) = “sin (konz cos0) g 
—ip sin (kọnz cos9) cos (konz cos@) 
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With, p = [Ecos 0 


The characteristic matrix for two layers is given as: 


M (z1) = Mo(zo)M1 (21 — Zo) 
M (z2) = M1 (z1)M2 (z2 — z1) 
Qo = M1 (21)Q (z1) (5) 


Where, Q (z1) = M2 (22 — z1)Q (22). 
The characteristic matrix of a jth layer is given by: 


1 ~~ kon;ôz; cos 8; 
M; = Bj AS (6) 
—ip;konjôzj cos 8; 1 
The characteristic matrix of the N layered media is given as: 
MS aed 1 ce 


With, A = et pjnj6z; cos 6; and B = yet T ôz; cos 6;. For an N layered media, this CMM is calculated as: 
j 


m m 
M= ie 7 | = (HL)" (8) 


The transmission (t) and reflection (r) coefficient is obtained as: 


t= Tac 2pı 
(mi1+mM12P1)p1+(M21+M22P1) 
pa R _ (mı1+mM12Pp1)p1-(M21+mM22P1) (9) 
A  (miıı+mı2p1)p1+(M21+M22P1) 


E E 
Here, pı = [= cos 6, and p; = [cos O, 


Where, €,, 44 and €,, 4, represents the electric permittivity and magnetic permeability of medium 1 and 
medium | respectively. The reflectivity (R) and transmissivity (T) are given as: 


R= Ir/? 
Tale? (10) 


4. PROPOSED SENSOR DESIGN 

The goal of the proposed work is to develop a theoretical design, simulation, and analysis of a 
distributed Bragg’s mirror for the detection of thyroid cancer. When the optical thickness of these two materials 
is equal to quarter wavelength of the incident electro-magnetic (EM) wave then the layered structure of two 
different materials exhibits a high reflectivity. A cavity can be created by removing any one of the layer (low or 
high layer). This makes the structure to allow one particular frequency of light to propagate through it. DBR 
has an alternate low and high RI materials, the optical thickness of each layer is quarter wavelength of 
incident light. The RI of L and H layer is represented as n, and ny respectively. The selected materials used 
for L and H layers are ZnS and Ge respectively. Silicon material is used as a base material over which this 
multilayer structure can be grown. Light incident on this layered structure undergoes multiple reflections in 
this media. The incident light is normal to the structure and is TE polarized. With the introduction of defect 
cavity in the middle layer of DBR structure, a sharp transmission of electromagnetic wave can be observed in 
the stopband region. Figure 3 shows the structure of defected multilayer cavity. A central cavity is 
sandwiched between N = 6 layer on either side. 
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Figure 3. Proposed fabry perot microcavity with two bragg’s mirror on either side with N = 3 


Above structure can be mathematically written as [27]. 


M, M, 


A lie Ma 


= (MM,)" Mp (MyM_)" (11) 


Where, M,, My and Mp represents L, H, and D layer the characteristic matrix respectively. 

The defect layer geometrical thickness is iterated carefully and finally chosen as 4 um. 
The transmission of resonant mode through the cavity is observed clearly for the normal and cancerous cells. 
In this work, the proposed micro-cavity structure is used to fill sample under test. The presence of thyroid 
cancerous cells can be detected based on the resonating wavelength of the light through the structure as soon 
as the sample is available at the central defect cavity. With the test sample in the central defect cavity, the 
incident light changes its resonating wavelength and there is a considerable wavelength shift for normal and 
cancerous cells. Hence, the presence of thyroid cancer cells can be detected in the given sample. In recent 
years, RI based detection is the widely used technique for the detection of cancerous cells. In this work, the 
RI of normal cells is taken as 0.8 and that of cancerous cells is taken as 0.9 [28]. Form this data it is clear that 
the RI of cancerous cells is higher than that of normal cells. 
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Figure 4. Normal cell transmission spectrum Figure 5. Thyroid cancer cell transmission spectrum 


As the incident light interacts with the sample in the central cavity and resonate with certain 
frequency. The presence and absence of thyroid cancer cells can be distinguished by observing this resonant 
frequency shift. The resonance spectrum for normal and cancerous cells is as shown in the Figure 4 and 
Figure 5 respectively. There is a significant shift between these two spectrum as shown in the figures. 

The resonant frequency of normal cell is 10.9042 um and that of cancerous cell is 11.187 um. This 
shift in the resonant wavelength calculated to determine the presence of cancerous cells in the given test 
sample. The wavelength shift is significant and is 0.2828 um. 
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In Table 1, we have shown the refractive index of normal and thyroid cancer cells [28]. The value of 
refractive index for normal cells is less than thyroid cancer cells. As in beginning of cancer, the deposition of 
protein particle in cells is starts so the value of refractive index of cells become larger. Also, the value of the 
refractive index of final stage cancer is larger than initial stage cancer. When the defect cavity is filled with the 
analyte containing cancerous and normal cells there is a significant shift in the resonating wavelength as shown in 
the combined Figure 6. Shift in the resonating wavelength for cancer and normal cells can be measured easily. 
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Figure 6. Combined transmission spectrum of normal and thyroid cancerous cells 


Table 1. Refractive index of normal and thyroid cancer cells [28] 
Cells under consideration Refractive index (RI) 
Normal 0.8 
Cancerous 0.9 


5. SENSOR PARAMETERS 

The performance and accuracy of the sensor are the important parameters of sensor for sensing 
applications. Here, we have calculated two major sensor parameters like sensitivity and quality factor [12], [13]. 
Both parameters such as sensitivity and quality factor are most important as sensor design point of view, so we 
calculate both parameters in this research work and tabulated in Table 2. The sensitivity of proposed biosensor 
is 2.828 (um/RIV) and Q factor is 3729. 


5.1. Sensitivity 

The sensing quality of a sensor is measured by using sensitivity of the sensor and is the ratio of 
resonant wavelength change to RI change of the sample. This is an important parameter of an optical sensor 
which describes the capability of the sensor to detect the desired component in a sample. 


Spy = E (12) 


5.2. Quality factor 

It is the ability of the sensor to detect malignant cells in the test sample. It is the ratio of resonant 
wavelength shift to the full width half maximum (FWHM) of a transmission spectrum. Q factor of an optical 
sensor describes the accurancy with which the desired component of an analyte can be measured. 


a 
= ores 13 
Q î1/2 ( ) 


Table 2. Proposed biosensor performance parameters 
RI of normal cells — Rl of thyroid cells Sensitivity (um/RIU) Q factor 
0.8 0.9 2.828 3729 


In order to justify the sensing ability of the proposed sensor, comparison table is shown in Table 3. 
It is clear from the table data that, the proposed sensor is having good sensing ability and Q factor as 
compared with the recently proposed works. So, in this way we can say that proposed sensor will give the 
better results than other sensors. 
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Table 3. Comparing proposed sensor with the recently reported works 


Serial number Sensor design Sensitivity (nm/RIU) Q factor Reference 
1 1D nano composite PhC [2019] 43 - [29] 
2 Ultra-sensitive PhC sensor [2019] 2200 10° [27] 
3 Square lattice defect-based PCW sensor [2020] 2360.12 99.765 [30] 
4 Proposed work 2828 3729 — 


6. CONCLUSION 

A theoretical design of one-dimensional Bragg’s reflector micro-cavity structure is proposed for the 
detection of thyroid cancer cells in the given test sample. For better sensor sensitivity an incident light is 
chosen in the mid-infrared range of frequency. On either side of the central defect cavity, 3 pairs of alternate 
low and high RI layers are present. Refractive index of normal and thyroid cancer cells is used to measure the 
sensor performance. The defect is filled with test sample and analysed the interaction of incident light with 
the sample. The resonant shift in wavelength for these cells is measured. The presence of thyroid cancerous 
cells is detected based on this wavelength shift. The sensor shows good Q factor and sensitivity of 3729 and 
2828 nm/RIU. Obtained results concludes that the proposed sensor can be suitable for medical applications 
for the quick and accurate detection of thyroid cancer cells. 
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